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We study the soft mode along the critical line in the phase diagram with the 
tricritical point, using the Nambu-Jona-Lasinio model. At the critical point with 
finite quark mass, the ordering density becomes a linear combination of the scalar, 
Q ' quark number and energy densities, and their susceptibilities diverge with the 

' same exponent. Based on the conservation law, it is argued that the divergent 

susuccptibility of a conserved density must be accompanied by a critically— slowing 
hydrodynamic mode. The shift of the soft mode from the sigma meson to the 
hydrodynamic mode occurs at the tricritical point on the critical line. 



1. Introduction 

The phases of QCD have been explored by investigating the behavior of 
the quark condensate, the Polyakov loop and the color superconducting 
gap as functions of the temperature (T) and the quark chemical potential 
(/i). Among various possibilities on the phase structure, existence of a 
critical point (CP) as an endpoint of the first-order phase boundary has 
been theoretically suggested and discussed in the literature 1,2,3 . 

When we extend the phase space by taking the mass (m) of the u and d 
quarks as the third variable, we can study this CP from the viewpoint of the 
phase space T-fi-m with a tricritical point (TCP). The static properties of 
this phase diagram with the TCP is well described by the Ginzburg-Landau 
(GL) effective potential of the quark condensate a expanded up to the a 6 
term. First in the case of exact chiral symmetry m = the T—\i plane 
must be devided into two domains of the symmetric and broken phases 
with a boundary line. Although the order of the singularity of this line is 



1 
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Figure 1. Phase diagram of the NJL model. The critical line is drawn in a dashed line. 
The first-order phase boundary forms a surface shown by hatch. The TCP is indicated 
by X. 



unknown in general, we assume the TCP where the order of the singularity 
shifts from the 2nd to the 1st order. Next, when the quark mass takes small 
but non-zero value, the 2nd order line disappears and the assumed TCP 
becomes a usual CP (see Fig. 1). In this consideration, the relation of this 
CP to the chiral symmetry is rather obscure. 

Characteristic time scale of the system response becomes infinitely large 
at a CP, which is known as critical slowing down. The mode whose typical 
frequency vanishes at the CP is called "soft mode." In the conventional 
theory increase of this time scale of a soft mode is related to the divergence 
of the susceptibility at the CP. The most familiar example will be the sigma 
meson or the radial fluctuation of the quark condensate in the chiral critical 
transition 5 . 

In this talk we discuss the soft mode along the line of the CP within 
effective approaches to QCD. 6 ' 7 We will point out that the soft mode as- 
sociated with the CP at finite m must have hydrodynamic character. This 
argument is based on the conservation of the baryon number and energy 
densities, and therefore is very general. Near the TCP the critical mode 
is different between the symmetric and broken phases. Although we re- 
strict ourselves to the result obtained in the Nambu-Jona-Lasinio (NJL) 
model here, we should stress that the same result can be reached within 
the time-dependent GL approach as well. 7 
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Figure 2. Effective potentials Q with two ordering densities, cr-p and a— s at critical 
poinsts, chiral CP (left), TCP (middle), CP (right). Those with the single ordering 
density (uppermost) are also shown. 



2. Effective potential 

We use the NJL model 8 ' 9 C = q(i@ — m)q + g[{qq) 2 + {qij5T a q) 2 ], in the 
mean field approximation ((qq) — cr=const, (qi-f5T a q) = 7r=0). The ther- 
modynamics is described by the effective potential, 

d 3 k 



/crh 
—-[E - Tln(l - n+) - Tln(l - n_)] 



4ff 



(23a) 2 , 



(1) 



where n± = {e^ E ^^ + E = VM 2 + k 2 , M = m - 2ga, and v = 

2NfN c = 2 • 2 • 3 = 12 with Nf and iV c the numbers of flavor and color, 
respectively. The true thermodynamic state is determined by the extremum 
condition, dVl/dcr = 0, and the corresponding grand potential is f2(T, \i, m). 
We define the model with the three-momentum cutoff A and with the 
coupling constant gA 2 = 2.5 which allows the TCP. In the following, all the 
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dimensionful quantities are expressed in the units of A. 

It is useful to difine the GL effective potential tl with two order pa- 
rameters in studying the true flat direction. 7 It is numerically constructed 
as shown in Fig. 2 at the chiral CP with (T,//)=(0.3419, 0.3), the TCP 
(0.20362, 0.49558) and a CP (0.1498, 0.5701) with m = 0.01 in the units 
of A. 6 The susceptibilities Xij — —yd 2 tt/didj = T,p,m) are equal 
to the inverse of the curvature matrix of the GL effective potential at the 
extremum point. Therefore the divergence of the susceptibilities at the crit- 
ical point is related to the appearance of a flat direction in the GL effective 
potential. 

At the chiral CP, the flat direction must be the scalar density a due 
to symmetry (Fig. 2). The susceptibility of a is divergent while those of 
the quark number density p and the entropy density s remain finite in the 
mean field approximation. 

The <7 2 and <r 4 terms in the potential (1) vanish at the TCP. This fact 
results in the large fluctuation along the potential valley of f2 with two 
ordering densities (Fig. 2). The fluctuations of p and s become divergent 
at the TCP approached from the broken phase. a If we define a potential 
of e.g., the quark density p by eliminating a with df2(cr, p)jda = 0, we see 
that the potential is flat on the lower density side of the critical density p t 
at the TCP. 7 

On the other hand, at the CP with the explicit breaking m ^ the 
proper ordering direction becomes a linear combination of a, p and s 
(Fig. 2). 10 ' 7 ' 11 The susceptibilities of these densities diverge with the same 
exponent since all of them involve a fraction of the critical fluctuation of 
the proper ordering density. Here the a direction is no longer special. One 
may choose equally well any of these densities as the ordering density in 
the static GL potential. 

3. Susceptibility and spectral density 

The susceptibility is obtained in the g— limit of the response function, which 
allows us to express the susceptibility as a sum of the spectral density: 

Xij = X*j(0,q^ 0) = hm / J - . (t,j = m,/i,T) (2) 

From this expression we see that divergence of the susceptibility is caused 
by spectral enhancement at ui = or mode softening, provided that the 



a Note that the fluctuation of s is a linear combination of those of the quark number 
density and the energy density. 
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spectral density 2Imxy (w,q) itself is integrable. 

At the chiral CP, the divergence of the scalar susceptibility is generated 
by softening of the sigma meson mode, which is the chiral partner of the 
pion mode. However, we should recognize that there is no symmetry reason 
to expect the massless sigma at the CP with explicit symmetry breaking 
due to to 7^ 0. 

There is a strong constraint on the spectrum in the fluctuations of the 
conserved quantity: the modes contributing to the susceptibility of a con- 
served quantity have to be hydrodynamic, that is, the typical frequency 
must vanish as q — > 0. Physically this is a consequence of the existence of 
the current j such that d t p + V • j — for e.g., the quark number density. 
Using the fact that the conserved density operator is commutative with the 
total Hamiltonian and the fluctuation-disspation theorem, we can express 
the susceptibility in another form, 8,7 

A 3 q-0 J 2tt 1 - e-^ V ' 

These two expressions (2) and (3) coincide with each other if and only if 
lim q ^o 2Imxij{oJ, q) = 2Tr5(Lo)uxij- Hence, when the susceptibility of a 
conserved quantity diverges, there must be a hydrodynamic mode which 
shows critical slowing. 

Finally we remark that the w-limit of the response function Xij( w ~~ * 
0, 0) has no contribution from the hydrodynamic mode spectrum. 



4. Soft modes in the NJL model 

The response functions in the random phase approximation are written as 

Xij(m,q) = 11^(^4, q) +U im (iq 4 ,q)- 1 -2gU mj (iq 4 , q). 

1 - lgYY mm (iq 4 , qj 

(4) 

Here polarization functions are defined with the imaginary-time quark 
propagator S{k) = l/(Jt + M) as 

n ij (%,q) = - / — -3T tr fcI) S(k)TS(k - q)T' , (5) 

n— — 00 

where q A = 2lirT (I e Z), k = (k, fc 4 + ifi) with fc 4 = -(2n + 1)ttT, T is 
an appropriate Dirac matrix, and the traces are taken over the flavor, color 
and Dirac indices. T = 1 for the scalar, ^74 for the baryon number, and 
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(a) (b) 

Figure 3. Spectral functions of the scalar channel above (Left, T/A = 0.35) and below 
(Right, T/A = 0.339) the chiral transition point with fi/A = 0.3 in the uj-q plane. 
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Figure 4. (a) Physical processes contributing to the spectrum with the detailed balance, 
(b) Analytic structure of the response function. The a meson and p-h poles (x) locate 
in the unphysical Riemann sheet. 



Hmf for (3 with 

Hmf = -i\r V +M + ifi l4 . (6) 

We deal with the response function of Hmf instead of the entropy because 
the entropy has no microscopic expression. The real-time response func- 
tion is obtained from the imaginary-time propagator through the usual 
replacement iq^ — > qo + it in the final expression. 

The collective modes arc generated by the infinite sum of the bubble dia- 
grams in the NJL model. The spectral function p mm (uj, q) = 2Imx mm (w, q) 
of the scalar response function 5 just above and below the chiral critical point 
with fixed fi = /i c are shown in Fig. 3. The spectrum in the time-like region 
comes from the quark pair creation/annihilation while the spectrum in the 
space-like region is due to the absorption/emission of the scalar fluctuation 
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Figure 5. Spectral functions of the scalar and the quark number susceptibilities at the 
CP with m/A = 0.01. 



by a quark or an anti-quark. The latter particle-hole (p-h) process gives 
rise to the Landau damping of collective motion in medium. These physical 
processes are schematically shown in Fig. 4 (a). In the unphyscal Riemann 
sheet of ui (Fig. 4 (b)), we found two kinds of complex poles corresponding 
to the collective excitations of these processes. We shall call here the pole 
related with the time-like spectrum the sigma meson and the pole for the 
space-like one the p-h mode. 

Approaching the chiral CP from the symmetric phase, we see in Fig. 3 
(a) that the sigma meson mode becomes soft. Meanwhile the p-h mode 
does not show any particular enhancement. From the broken phase, on the 
other hand, the mass gap of the sigma meson mode is vanishing and the 
p-h mode spectrum also gets stronger at q = (Fig. 3 (b)). 

At the critical point with to/ A = 0.01 both the scalar and quark-number 
susceptibilities diverge. The scalar spectral function is shown in Fig. 5 (a). 
In this case the sigma meson has finite mass gap, which is set by the quark 
mass as ~ 2M ~ to 1 / 5 . The p-h mode has the hydrodynamic character 
(cj — > as q — > 0). We note that the strength of this mode is strongly 
enhanced at this critical point. In Fig. 5 (b), we show the specral function 
of the quark-number response function. It is clear that the divergence is 
caused solely by the p-h mode spectrum with hydrodynamic character, 
which is consistent with the general argument given in Sec. 3. 
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Figure 6. Ratio (7) of the spectral contribution along the chiral critical line in Fig. 1 
approached from the broken phase. R — > 1 toward the TCP. 



5. Change of the spectral contribution along the critial line 

We know that the sigma meson mode is the soft mode associated with the 
chiral CP while we have seen that the p-h mode shows the critical slowing 
at the CP with finite quark mass m. Let us study the change of the critical 
eigenmode along the critical line as shown in a dashed line in Fig. 1, with 
defining the relative weight of the spectral contributions of the p-h mode 
to the total spectrum by 



Here we used the fact that the difference between the q- and w-limits of 
Xmm(w, q) stems from the spectral contribution of the hydrodynamic mode. 
From the argument in Sec. 3, this ratio must be unity for the susceptibility 
of a conserved quantity, which can be confirmed explicitly. 

We show the numerical result of R in Fig. 6. The ratio vanishes if the 
chiral CP or the TCP is approaced from the symmetric phase, which means 
that the critical divergence is generated by the sigma mode without any 
contribution from the p-h mode. From the broken phase, on the other 
hand, we see the finite portion of the divergence comes from the p-h mode 
spectrum in the chiral critical transition — even in the fi = case. Ap- 
proaching the TCP along the critical line, we see that the contribution from 
the p-h mode increases and eventually saturates the spectral sum (2). 

First one may ask why the p-h mode gives no contibution in the sym- 



R = 



Xm m (0,0+) -y mm (0 + ,0) 

Xmm(0,O+) 



(7) 
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metric phase. This is because that the chirality and hclicity is the same for 
a massless quark and that the chirality flip in the scalar coupling requires 
the finite momentum transfer q. In the broken phase, a massive quark with 
definite chirality has both helicity components, and therefore the p-h mode 
spectrum remains contributing to the susceptibility in the q — > limit. 

Next we note that the mixing of this p-h spectrum in the susceptibility 
of p and s is possible in the broken phase through the coupling proportional 
to the condensate a. This mixing gives rise to the finite gap of these 
susceptibilities across the chiral critical line. This fact implies that the p-h 
mode contribution must be of order of 1/c 2 — > oo in the scalar channel, to 
cancel the a 2 factor from the coupling. 

At the TCP approached from the broken phase with fixed /i, \ mm di- 
verges as 1 /a 4 while x w an d Xtt blow up as 1/cr 2 , which is easily confirmed 
within the Ginzburg-Landau approach. Only the p-h mode spectrum with 
hydrodynamic character can cause this divergence in the NJL model. The 
sigma meson mode in the NJL model provides the singularity of order 1 jo 2 
to Xmm at the TCP approached from the broken phase. From the symmeric 
phase, where a = 0, there is no contribution from the p-h mode and the 
critical divergence at the TCP is completely provided by softening of the 
sigma meson. The x w and xtt are finite there. 

Along the critical line with finite quark masses m ^ 0, x w and xtt 
as well as Xnn diverge with the same exponent. Generally these divergence 
must come from softening of a hydrodynamic mode in the system because 
the baryon number and the energy are conserved quantities. We can prove 
this within the NJL model as well as the time-dependent GL approach. 

6. Summary 

We have seen that, unlike at the chiral CP, the ordering density at the CP 
with finite quark mass m is a linear combination of the scalar density, the 
baryon number density and the energy density, and that the susceptibilities 
of these density diverge there. Since the susceptibility of a conserved density 
solely comes from the hydrodynamic spectrum, the associated critical soft 
mode must be hydrodynamic. We identified the p-h mode in the NJL model 
as this critical mode. Recently it is explicitly argued that the dynamic 
unversality class of this point is the same as the liquid-gas critical point 11,12 . 

Experimentally, it is worthwhile to study the fluctuation of the con- 
served densities. 1 The correct evolution equation for the correlation length 
must be hydrodynamic one, which is slower than the sigma like motion. 
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In this sense the growth of the correlation length in the heavy ion events 
passing by the CP acquires a renewed interest. 13 

Along the critical line in the phase space of T-fi-m, we have studied the 
changeover of the associated soft mode from the sigma meson mode at the 
chiral CP to the p-h mode at the CP with m ^ 0. This shift occurs at the 
TCP, where the critical soft mode is different between the symmetric and 
broken phases. The dynamic classification of this TCP will be theoretically 
interesting. 7 ' 11,12 

There is a speculated phase diagram of QCD in the space of T—^i-m— 
771,5, 14 in which we see several critical lines and surfaces. It seems important 
to keep the hydrodynamic viewpoint in mind when we study these critical- 
ities. 
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